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ABSTRACT 

The Cygnus OB2 Association is one of the nearest and largest collections of massive stars in the 
Galaxy. Situated at the heart of the “Gygnus X” complex of star-forming regions and molecular 
clouds, its distance has proven elusive owing to the ambiguous nature of kinematic distances along 
this ^ ~ 80° sightline and the heavy, patchy extinction. In an effort to refine the three-dimensional 
geometry of key Cygnus X constituents, we have measured distances to four eclipsing double-lined OB- 
type spectroscopic binaries that are probable members of Cyg OB2. We find distances of 1.33 ± 0.17, 

1.32 ± 0.07, 1.44 ± 0.18, and 1.32 ± 0.13 kpc toward MT91 372, MT91 696, CPR2002 A36, and 
Schulte 3 respectively. We adopt a weighted average distance of 1.33±0.06 kpc. This agrees well 
with spectrophotometric estimates for the Association as a whole and with parallax measurements of 
protostellar masers in the surrounding interstellar clouds, thereby linking the ongoing star formation 
in these clouds with Cyg OB2. We also identify Schulte 3C (09.5V), a 4" visual companion to the 
4.75 day binary Schulte 3(A-|-B), as a previously unrecognized Association member. 

Subject headings: Techniques: radial velocities — (Stars:) binaries: general — (Stars:) binaries: 

spectroscopic — (Stars:) binaries: {including multiple) close — Stars: early-type 
— Stars: kinematics and dynamics 


1. INTRODUCTION 

The Cygnus OB2 Association (also known as VI Cygni) 
harbors the largest number of O stars among nine OB 
associations within the “Cygnus X” (Piddington & Min- 
nett 1952) complex, and it is one of the richest OB as¬ 
sociations in the Galaxy. Cygnus X is a conspicuously 
bright region in radio. X-ray, and infrared maps of the 
Galactic Plane, covering nearly ten degrees near lon¬ 
gitude £=80° (review by Odenwald & Schwartz 1993). 
Discovered as an overdensity of bright, reddened, early- 
type stars (Miinch & Morgan 1953; Johnson & Mor¬ 
gan 1954; Schulte 1956, 1958), Gyg OB2 boasts ^60 O 
stars (Wright et al. 2015), perhaps as many as ~2600 B 
stars (Knodlseder 2000), and some of the most massive 
stars known (Massey & Thompson 1991; Gomeron et al. 
2002). These are all veiled by heavy and variable inter¬ 
stellar extinction (Av= 4-12 mag; Torres-Dodgen et 
al. 1991; Massey & Thompson 1991; Hanson 2003). In 
recent years, large surveys have produced an unprece¬ 
dented panchromatic view of the interstellar environ¬ 
ment that produced Gyg OB2. These include a five-band 
mid-infrared survey using the Spitzer Space Telescope 
(Beerer et al. 2010), maps of the molecular clouds in 
several GO transitions (Schneider et al. 2006; Leung & 
Thaddeus 1992), and maps of the dust continuum at 
1.2 mm (Motte et al. 2007). O stars in Cyg OB2 were 
the first known stellar sources of X-ray emission (Barn- 
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den et al. 1979; Waldron et al. 1998), and several of 
the luminous members of the Association are prototyp¬ 
ical examples for whole classes of massive star phenom¬ 
ena, such as the extremely luminous early-B supergiant 
Schulte 12 (MT91 314 in the numeration of Massey & 
Thompson 1991), the non-thermal radio source Schulte 
9 (MT91 431; Abbott et al. 1984; Naze et al. 2012), 
and the 03If* star Schulte 7 (MT457; Walborn 1973). 
As one of the nearest regions of massive star formation, 
Cyg OB2 also permits detailed investigations of OB stars 
and their interplay with the interstellar medium, such as 
the stellar bowshocks associated with high-velocity mem¬ 
bers (Comeron & Pasquali 2007; Gvaramadze 2007; Kob- 
ulnicky et al. 2010). 

Despite its role as a nearby laboratory of massive star 
astrophysics, the distance to Cyg OB2 has proven elusive. 
In part, this is the result of variable extinction, which 
complicates spectrophotometric distances. Additionally, 
its location is at £ = 80°, where the relation between 
radial velocity and distance is not only double-valued 
but poorly-defined owing to velocity crowding along the 
tangent point of the local spiral arm (e.g., see Dame & 
Thaddeus 1985; Dame et al. 2001). Spectrophotomet¬ 
ric distance measurements span a range that includes 
1.5 kpc (Johnson & Morgan 1954), 2.1 kpc (Reddish et 
al. 1966), 1.7 kpc (Massey & Thompson 1991), 1.7 kpc 
(Torres-Dodgen et al. 1991), and 1.45 kpc (Hanson 2003). 
Additionally, studies of the eclipsing, double-lined spec¬ 
troscopic binaries Cyg OB2-B17 (Stroud et al. 2010) and 
V382 Cyg (Ya§arsoy & Yakut 2012) give distances of 1.5 - 
1.8 kpc and 1.466 ± 0.076 kpc respectively. This range 
places Cyg OB2 within the local spiral arm or spur but 
is insufficiently precise for measuring key physical pa¬ 
rameters such as luminosities and mass loss rates. Rygl 
et al. (2012) used radio-wave very-long-baseline interfer- 
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ometry (VLBI) to measure parallaxes of five pre-main- 
sequence water and OH masers within Cygnus X and 
found that four of the five had distances between 1.3 and 
1.5 kpc with typical uncertainties of 0.1 kpc. Zhang et 
al. (2012) also used VLBI techniques to measure a paral¬ 
lax distance of 1.61±0.12 kpc to the red hypergiant NML 
Cygni, and they concluded that it may lie on the far side 
of the Cyg OB2 association. Dzib et al. (2012) reported 
radio VLBA observations giving a parallax distance of 
l-65lg;44 kpc to the colliding-wind binary Schulte 5, and 
they infer a probable distance of 1.3 - 1.4 kpc based on 
estimates of the orbital parameters in this quadruple sys¬ 
tem. While these parallax determinations are the most 
precise distances yet published, Cygnus X probably has 
a signihcant spread along the line of sight. Because of 
this, individual OB associations and star forming regions 
may lie anywhere within this depth. Thus, independent 
distance measurements of individual features are needed 
to help refine the line of sight structure of this complex, 
in particular for Cyg OB2. 

In this contribution, we use four eclipsing double-lined 
spectroscopic binaries (SB2s) investigated as part of the 
Cygnus 0132 Radial Velocity Survey to calculate an in¬ 
dependent distance to the Association. Papers I - VII in 
this series (Kiminki et al. 2007, 2008, 2009, 2012; Kiminki 

6 Kobulnicky 2012; Kobulnicky et al. 2012, 2014) de¬ 
scribe prior results of the Cygnus OB2 Radial Velocity 
Survey which measured the massive binary characteris¬ 
tics (i.e., binary fraction, distribution of periods, mass 
ratios, eccentricities) for more than 100 massive stars in 
a single cluster/association having a common formation 
environment. Paper V (Kiminki & Kobulnicky 2012), 
in particular, uses these basic data to infer the intrinsic 
distributions of massive binaries, concluding that: the 
fraction of massive stars having companions may be as 
high as 90%, 45% of these multiples are likely to interact 
at some point, there exists an excess of short-period 4 - 

7 day systems relative to 7 - 14 day systems, and that 
unresolved secondaries contribute ~16% of the light in 
young stellar populations. Paper VII uses the extensive 
repository of orbital parameter information for the 48 
known massive binaries to show that no single power law 
provides a statistically compelling prescription of the cu¬ 
mulative orbital period distribution, and that a flattening 
of the distribution at P > 45 d indicates either a lower 
binary fraction or a shift toward low-mass companions 
among long-period systems. 

Our present dataset includes optical spectra and 
broadband photometry of four eclipsing double-lined 
spectroscopic binaries: MT91 372 and MT9I 696 
(nomenclature of Massey & Thompson 1991), 
CPR2002 A36 (nomenclature of Comeron et al. 
2002, also known as RLP357 in Reddish et al. 1966), 
and Schulte 3 (nomenclature of Schulte 1958, also 
known as RLP920 in Reddish et al. 1966). Addition¬ 
ally, we make use of the Northern Sky Variability 
Survey (NSVS) photometry (Wozniak et al. 2004) for 
CPR2002 A36 and Schulte 3. In Paper III (Kiminki et 
al. 2009), we presented a single-lined orbital solution 
and NSVS photometry with a period of 2.22 days for 
MT91 372 (B2?-|-early B). Rios & DeGioia-Eastwood 
(2004) reported the SB2 (09.5V-|-early B) nature of 
MT91 696 and estimated a photometric period of 1.46 


d. A rehned photometric period of 1.46919 ± 0.00006 d 
appears in Souza et al. (2014). We presented the 
first spectroscopic orbital solution in Paper IV with 
a similar period (Kiminki et al. 2012). Otero (2008a) 
identify the eclipsing nature of CPR2002 A36 using 
NSVS photometry. Its SB2 nature was noticed by 
Hanson (2003), and we presented the orbital solution 
for this P=4.67 d system in Paper HI (Kiminki et al. 
2009). Schulte 3 was discovered to be a double-lined 
(06IV-I-09HI) spectroscopic binary with P=4.74 d in 
Paper H of this series (Kiminki et al. 2008). In this 
work, we analyze the light curve of this eclipsing system 
for the first time. 

Section 2 of this work describes the new photometric 
and spectroscopic data obtained for the purpose of cre¬ 
ating joint velocity curves and light curves for these four 
eclipsing double-lined systems. Section 3 reports the dis¬ 
covery of additional visual companions at small angular 
separations from Schulte 3 and MT91 696; knowledge 
of these third-light contributions is essential for accurate 
luminosities and distances. Section 4 provides details of 
the combined light-curve/velocity-curve modeling used 
to infer the distance to each system. Section 5 summa¬ 
rizes the new eclipsing binary distance to Cyg OB2 in 
relation to measurements of other targets in this region. 
All radial velocities reported here are in the Heliocentric 
reference frame. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. Spectroscopic Data 

Optical spectroscopic observations of the four targets 
were obtained as part of the Cygnus OB2 Radial Velocity 
Survey on numerous dates between 2001 and 2011 using 
a variety of telescopes and instruments as described in 
Papers I - IV of this series. Several additional obser¬ 
vations were obtained in 2012 October (on MT91 696 
and Schulte 3) and 2013 May - June (on MT91 372) us¬ 
ing the Wyoming Infrared Observatory (WIRO) Longslit 
spectrograph with an e2V 2048^ CCD as the detector. 
A 2000 1 mm“^ grating in first order yielded a spec¬ 
tral resolution of 1.25 A near 5800 A with a 172x100" 
slit. The spectral coverage was 5250 - 6750 A. Expo¬ 
sure times ranged between 1200 s and 2 hours in mul¬ 
tiples of 600 s depending on target brightness, current 
seeing (172 - 3" FWHM) and cloud conditions. Re¬ 
ductions followed standard longslit techniques, including 
flat fielding from dome quartz lamp exposures. Copper- 
argon arc lamp exposures were taken after each star ex¬ 
posure to wavelength calibrate the spectra to an rms 
of 0.03 A (1.5 km s“^ at 5800 A). Multiple exposures 
were combined yielding final signal-to-noise ratios (SNR) 
typically in excess of 100:1 near 5800 A. Einal spec¬ 
tra were Doppler corrected to the Heliocentric velocity 
frame. Each spectrum was then shifted by a small addi¬ 
tional amount in velocity so that the Na I D AA5890,5996 
lines were registered with the mean Na I line wavelength 
across the ensemble of observations. This zero-point cor¬ 
rection to each observation is needed to account for ef¬ 
fects of image wander in the dispersion direction when 
the FWHM of the stellar point spread function was less 
than the slit width. Because of these inevitable slit- 
placement effects on the resulting wavelength solutions, 
radial velocity (RV) standards were not observed. Rel- 
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ative velocity shifts were generally less than 6 km s“^, 
comparable to the magnitude of the random velocity un¬ 
certainties. 

We measured an initial (final for CPR2002 A36) ra¬ 
dial velocity for each spectrum obtained at WIRO us¬ 
ing Gaussian fits to the He I A5876 line via the same 
method outlined in Kobulnicky et al. (2012), adopting 
a rest wavelength of 5875.69 A measured in model stel¬ 
lar spectra with static atmospheres (TLUSTY; Lanz & 
Hubeny 2003; Hubeny & Lanz 1995). Our fitting code^ 
fixes the Gaussian width and depth at the mean value 
determined from all the spectra (after rejecting outliers), 
and it solves for the best-fit line center. In the case of an 
SB2, the code fits a double-Gaussian profile where the 
widths and depths have been fixed independently using 
observations obtained near quadrature orbital phases. 

In this work, we used the initial radial velocities to dis¬ 
entangle the component spectra of MT91 372, MT91 696, 
and Schulte 3 using the method of Gonzalez & Levato 
(2006).® One of the strengths of this method is that the 
radial velocities can be refined via cross correlation after 
each iteration (i.e, cross-correlating the residual spectra 
with the resultant component spectrum as the template). 
The cross-correlated velocities generally have smaller un¬ 
certainties and utilize more lines, leading to a more pre¬ 
cise measurement of the stellar velocity and the com¬ 
puted binary systemic velocity. Systemic velocities that 
are based solely on one hydrogen or helium line, such as 
He I A5876 with CPR2002 A36, may be blueshifted rela¬ 
tive to the true systemic velocity owing to velocity con¬ 
tributions from strong winds. These contributions are 
strongest in very early-0 stars and evolved stars. How¬ 
ever, radial velocities for Schulte 3 are measured from 
many lines in the blue and red portions of the spec¬ 
trum. The radial velocities of MT91 372 and MT91 696 
are based on several hydrogen and helium lines between 
5400 A and 6700 A and the stars are neither early-0 or 
evolved. 

Because of the more sophisticated approach to mea¬ 
suring the radial velocities in this work, the newly de¬ 
termined velocities will vary slightly from our previous 
works. This will be most evident for Schulte 3 and 
CPR2002 A36, which now have smaller 0-G velocity 
residuals, indicating that a superior orbital solution has 
been obtained with the use of new data and methods. 
Table A.l in the Appendix lists the Heliocentric Julian 
dates (column 1), orbital phase (column 2), radial veloc¬ 
ities and uncertainties (columns 3 and 5), and the 0-C 
residuals (columns 4 and 6) of each component for each 
of the four eclipsing SB2 systems. 

2.2. Photometric Data 

Photometry of MT91 696 was obtained on nine nights 
between 2012 October 17 and November 10 at the Uni¬ 
versity of Wyoming 0.6 meter Red Buttes Observa¬ 
tory (RBO) using a 1024x1024 Apogee Alta U47 GGD 
through a V filter. Data on Schulte 3 was obtained on 
14 nights at RBO between 2012 October 17 and 2012 
December 21, and data on GPR2002 A36 was obtained 

® We use the robust curve fitting algorithm MPFIT as imple¬ 
mented in IDL (Markwardt 2009). 

® Disentangling is not performed for CPR2002 A36 on account 
of having an insufficient number of high-quality spectra. 


on four nights between 2012 December 10 and 2012 De¬ 
cember 22. The 0'.'54 pixel“^ scale at RBO yielded a 
9.0' field of view centered on each object. Photometry 
of MT91 372 was also obtained at RBO on seven nights 
between 2013 June 14 and 2013 July 8 using an Alta U16 
4096x4096 GGD binned 2x2, yielding 0'.'73 pixels over a 
24' field of view. The observations span 10 orbital peri¬ 
ods of MT91 372, 16 orbital periods of MT91 696, two 
orbital phases of CPR2002 A36, and 13 orbital periods 
of Schulte 3. Seeing varied between 3" and 5" FWHM. 
Sky conditions were predominantly non-photometric. Se¬ 
quences of 180 s, 120 s, and 30 s exposures were obtained 
for all systems, for several hours each night. Photo¬ 
metric conditions on 2012 October 18 (local) allowed us 
to obtain calibration images of several Landolt standard 
fields and fields in Gyg OB2 having UBV photometry in 
Massey & Thompson 1991, (MT91). We elected to per¬ 
form differential photometry using OB stars measured by 
Massey & Thompson (1991) in the U-band because of the 
lack of suitably red {B — V > 1.2) Landolt standards re¬ 
quired to match the colors of the heavily reddened OB 
star targets. Additionally, we observed each target in 
the UBV filters on the (mostly) photometric night of 
2012 December 09, interleaving target fields with nearby 
Cyg OB2 fields for purposes of tying our photometry to 
the MT91 UBV measurements. 

RBO images were reduced using standard procedures 
which included removing the zero level with bias expo¬ 
sures, removing dark current using a median of at least 
seven dark exposures scaled to the appropriate target 
exposure time, and dividing by a normalized flat field 
constructed from the median of at least seven twilight 
sky exposures having a minimum of 10,000 electrons 
per image. We performed aperture photometry using 
IRAF/PHOT with an aperture of 14 pixels (8") and a 
sky annulus extending from 14 to 19 pixels, in order to 
extract > 99% of the stellar flux for the target and 6-9 
reference stars in each field. The magnitudes of the tar¬ 
get stars were compared to the average of the compar¬ 
ison stars in each field to correct for (sometimes large) 
variations in atmospheric transmission. Differential pho¬ 
tometry of MT91 372 used nine nearby field stars for ref¬ 
erence, including four early type Cyg OB2 stars having 
similar color. The zero point of the photometric cali¬ 
bration relies upon stars from MT91 in the same held. 
The magnitudes of each comparison star were examined 
relative to the mean for the remainder of the stars, and 
in this way, we judge that the comparison stars were 
non-variable to the level of 0.015 mag over the course of 
our program. Photometric errors of the target star lie 
in the range 0.001 - 0.005 mag when considering only 
photon statistics as the dominant source of uncertainty. 
The exposure-to-exposure rms suggests that this is an 
appropriate estimate of the relative photometric uncer¬ 
tainty on any given night once brief time periods of heavy 
cloud cover were excised from the data. However, the rms 
difference of the comparison star mean magnitudes over 
all observing nights suggests that the night-to-night zero 
point is uncertain at the 0.010 mag level, which is small 
compared to the depth of eclipses over the orbital cycle. 

For Schulte 3 and GPR2002 A36, we included the 
NSVS photometry (Wozniak et al. 2004) in addition to 
the RBO photometry. The NSVS photometry were taken 
with four Apogee AP-10 imaging cameras with an until- 
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tered optical response of ~450 - 1000 nm and an effective 
wavelength close to Johnson R-band. The data were col¬ 
lected over the course of one year, from 1999 April 1 
through 2000 March 30 and exhibit an average photo¬ 
metric uncertainty of 0.01 mag. The NSVS data used 
here are converted from Modified-Julian Date (MJD) to 
Julian Date (JD). 

Tables A.2 - A.5 in the Appendix list the Helio¬ 
centric Julian date and the H-band magnitudes of the 
MT91 372, MT91 696, CPR2002 A36, and Schulte 3 sys¬ 
tems. Column 1 gives the Heliocentric Julian dates mi¬ 
nus 2,400,000, while columns 2 and 3 list the H-band 
magnitudes and Icr uncertainties calculated from pho¬ 
ton statistics only. Additional uncertainties at the level 
of 0.010 mag are present as a result of uncorrected at¬ 
mospheric transparency variations over the field of view 
and night-to-night zero point uncertainties. Only the 
first ten table lines are given to provide context and for¬ 
mat. The full table appears in the electronic edition 
as a machine-readable table. The Schulte 3 photome¬ 
try includes the contribution from the blended tertiary 
Schulte 3C discussed in the next section. Additionally, 
the data have been averaged into 2-minute intervals from 
the original 30-second exposure times. MT91 696 pho¬ 
tometry also includes a contribution from the blended 
tertiary MT91 696C, discussed in the next section. 

3. VISUAL COMPANIONS TO SCHULTE 3 AND 
MT91 696 

3.1. Schulte 3C 

Photometry of the double-lined eclipsing binary 
Schulte 3 was complicated by the presence of a nearby 
companion at 379 separation at position angle 215° from 
the brightest visual component. Schulte 3 (comprised of 
4.74 d eclipsing components “A” and “B”) and the com¬ 
panion, Schulte 3C, are blended in all of our RBO images, 
so our photometry using large apertures encompasses 
all stars. The companion is readily visible in UKIDSS 
Galactic Plane Survey images (Lucas et al. 2008) at 
R.A.=20:31:37.31 Declination=-t41:13:17. (J2000), but 
both Schulte 3(A-|-B) and Schulte 3C are saturated. The 
companion is not present in the 2MASS Point Source 
Catalog because of blending. It is not visible in the 
Digitized Sky Survey for the same reason. We there¬ 
fore conducted our own U-band imaging using the CCD 
guide camera on the WIRO-Longslit spectrograph (07147 
pixel”on the night of 2012 November 06. Several 3 
s exposures were obtained in 174 seeing within a few 
minutes of JD=2456238.55 (2012 Nov 7, 01:19 UT). Fig¬ 
ure 1 shows a greyscale representation of our U-band im¬ 
ages, illustrating the relative locations of Schulte 3(A-|-B) 
and Schulte 3C. Photometry using these images indicates 
a magnitude difference Amy=2.14 mag, meaning that 
Schulte 3(A-|-B) contribute 88% of the flux at V band, 
compared to 12% for Schulte 3C. Given a mean mag¬ 
nitude of V=10.30±0.02 on this date, derived from the 
phased light curve of Schulte 3 (§4.5), the companion has 
V~12.42±0.03. It also means the thirdlight-adjusted V- 
band maximum (between primary and secondary eclipse) 
for the Schulte 3(A-|-B) system, as derived from the 
phased light curve, is 10.31 mag. 

We obtained spectra of Schulte 3(A-|-B) and Schulte 3C 
using the WIRO Longslit spectrograph on the night of 
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Fig. 1.— Greyscale image of Schulte 3 (the eclipsing double-lined 
components a-|-b) and its ~09.5V companion, Schulte 3C. 


2012 October 21 in 2" seeing. Figure 2 plots the spec¬ 
tra of both components and labels some key spectral 
features of hot stars. The spectrum of Schulte 3C ex¬ 
hibits both He II A5411 and He I A5876 in the ratio 
EW 5411 /EW 5876 0.4, consistent with a temperature 

class near 09.5 (Kobulnicky et al. 2012). Ha appears in 
absorption. The combination of the spectral characteris¬ 
tics and the U-band magnitude are both consistent with 
an 09.5V at the distance of Cyg OB2. With the current 
data it is not possible to say whether the companion is 
physically associated with Schulte 3, although the strong 
similarity of the many diffuse interstellar features in the 
two spectra make it likely that they both lie at a similar 
distance. If associated, Schulte 3 is a triple stellar system 
with Schulte 3C at a projected separation of about 5200 
AU. 



5400 5600 5800 6000 6200 6400 6600 

Wavelength ( A ) 


Fig. 2.— Spectra of the 06V-h09III double-lined system 
Schulte 3(A-hB) (blue) and its 4'^ 09.5V companion, Schulte 3C 
(red). Labels mark key stellar features. Other features are inter¬ 
stellar. The excellent agreement between the interstellar features 
of both stars are consistent with a common distance. The two 
stars differ most notably in their He I A5876 and Ha line depths, 
reflecting the higher effective temperature of Schulte 3(A-|-B). 
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In the analysis that follows, we assume, but do not 
have the data to demonstrate, that Schulte 3C is non¬ 
variable, both photometrically and spectroscopically. We 
are confident that it cannot be the source of the photo¬ 
metric variability, detailed in the next section because 
the spectroscopic period of Schulte 3(A-|-B), as observed 
in isolation, matches the photometric period observed for 
the whole system. However, the present data do not rule 
out the possibility that Schulte 3C may be a multiple 
subsystem or may be photometrically variable. 

3.2. MT91 696C and 696D 

While no similarly bright companions are seen in 
UKIDSS images for MT91 696, measurements with 
the Hubble Space Telescope Fine Guidance Sensors 
(FGS) (Caballero-Nieves et al. 2014) reveal two astro- 
metric companions at 0'.'023 (MT91 696G) and Cf.'Sd 
(MT91 696D) separations. A Cyg OB2 distance of 
1.33 kpc implies minimum projected physical separa¬ 
tions of 31 AU and 1117 AU respectively, indicating that 
this is likely a triple (possibly quadruple) system com¬ 
posed of a spatially unresolved 1.4 day binary (compo¬ 
nents A & B) and two additional components (C & D). 
With Am=0.94±0.40 mag in the HST F583W band¬ 
pass, MT91 696C is the dominant source of excess light 
(MT91 696D is several magnitudes fainter), so we con¬ 
sider only this component in subsequent analysis of the 
light curve. The magnitude of MT91 696C would make it 
a probable mid-B star. As with Schulte 3, we assume, but 
do not have the data to demonstrate, that MT91 696C 
is non-variable, both photometrically and spectroscopi¬ 
cally. It is also not the source of photometric variability 
for the same reasons as Schulte 3C. Lastly, while we con¬ 
sider photometric contamination from MT91 696C, we 
assume negligible radial velocity contamination given a 
computed semi-amplitude of ~5.5 km s“^and total flux 
contribution of ^14% (see Section 4.4). 

4. ANALYSIS OF ECLIPSING BINARIES 

4.1. Modeling the Joint Light Curves and Velocity 
Curves 

In general, effective temperatures for the primary star, 
Ti, were held as fixed parameters based on the ratio of 
equivalent widths He II A5411/He I A5876 in our spectra. 
Figure 3 plots this EW ratio versus temperature as mea¬ 
sured in Tlusty NLTE model atmosphere spectra (Lanz 
& Hubeny 2003) and CMFGEN model spectra (Hillier & 
Miller 1998)’’ for four values of log(g), 4.0 - 3.25 (appro¬ 
priate for luminosity classes V - I). Between r=30,000 K 
and ~35,500 K there is a nearly linear relationship be¬ 
tween EW ratio and temperature. The solid line is a fit to 
the log(5')=4.0 - 3.75 models over the range T=27,500 - 
40,000 K, and the dashed line is a fit to the models with 
log(c,)=3.25 - 3.50 over the range r=27,500 - 37,500 K. 
Above these maximum temperatures, the linear relation 
becomes a poor approximation. Labels near the top of 
the panel give the corresponding spectral type according 
to the observational calibration of Martins et al. (2005) 
for dwarf and supergiant luminosity classifications. 

^ As measured from the grid of O star stellar mod¬ 
els dated 2009 June 24 from the webpage of J. Hillier, 
http://kookaburra. phyast.pitt.edu/hillier/web/CMFGEN.htm. 



Fig. 3.— Relation between equivalent width ratio He II 
A5411/He I A5876 versus temperature based on Tlusty Lanz & 
Hubeny (2003) and CMFGEN (Hillier & Miller 1998) NLTE model 
atmospheres for values of log( 3 ) appropriate to main-sequence 
(3.75 — 4.0) and giant/supergiant (3.25 - 3.5) late-O and early B 
stars. The expressions within the panel quantify this quasi-linear 
relationship over the range 27,500 K - 37,500 K for luminosity class 
V stars (solid line) and 27,500 K - 35,000 K for luminosity class 
I/HI stars (dashed line), where is temperature in 10® K. 

To obtain the most probable orbital period, we exam¬ 
ined each component’s CLEANed radial velocity power 
spectrum using an IDL® program written by A. W. 
Fullerton, which makes use of the discrete Fourier trans¬ 
form and CLEAN deconvolution algorithm of Roberts 
et al. (1987). Sufficient radial velocity and light curve 
data exist for each system so that periods are secure and 
free of aliases. Spectroscopic solutions were then com¬ 
puted with the Binary Star Combined Solution Package, 
BSCSP (Gudehus 2001). The spectroscopic orbital param¬ 
eters determined from BSCSP were used as initial param¬ 
eters within the eclipsing binary package, PHOEBE (Prsa 
& Zwitter 2005), which is based on the code of Wilson 
& Devinney (1971), to model the joint light curves and 
radial velocity curves. For the light curve parameters, 
we chose albedos and gravity darkening coefficients ap¬ 
propriate for radiative atmospheres (1.0) and limb dark¬ 
ening coefficients interpolated from van Hamme (1993). 
Gravity darkening and albedo were treated as fixed pa¬ 
rameters. We allowed each system one reflection from 
the opposite star and used the Kurucz atmospheric mod¬ 
els as the description of the star’s wavelength-dependent 
emergent flux (Kurucz 1979). 

We selected the appropriate PHOEBE geometrical con¬ 
figuration for each system based on the shape of 
the light curve (i.e., double-contact for GPR2002 A36 
and Schulte 3, and detached for both MT91 372 and 
MT91 696). Including only the spectroscopic data for 
each system, we used PHOEBE’s differential correction 
method to refine the spectroscopic solution while allow¬ 
ing period, eccentricity, separation, epoch of periastron, 
angle of periastron, systemic velocity, and mass ratio 
to vary. From this point, for the detached systems, we 
determined the roughly optimized photometric solution 
by varying inclination, surface potential, luminosity, and 

® The Interactive Data Language (IDL) software is provided by 
Exelis Visual Information Solutions. 
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secondary temperature, using only the photometric data. 
Schulte 3 and CPR2002 A36 show evidence of hot spots 
on the secondary component via the O’Connell effect in 
their light curves®, Struve-Sahade effect^® in their He I 
lines, and complex varying Ha emission. For the model 
of these two systems, we added one spot to the secondary 
and manually adjusted the spot temperature, size, and 
location in an iterative process with the differential cor¬ 
rections to obtain the roughly optimized solution. De¬ 
spite the evidence pointing to the spot being associated 
with the secondary (discussed in Sections 4.5 and 4.6), 
we also attempted placing the spot on the primary. For 
both systems, the spot on the secondary provided the 
better fit to the data. The spot parameters were ad¬ 
justed manually because of the degenerate nature of the 
solutions (i.e., spot size, spot temperature, and stellar 
effective temperature are degenerate when the analysis 
is limited to 1 - 2 photometric bandpasses). 

Once refined solutions for both the photometric and 
spectroscopic data were obtained separately, we opti¬ 
mized the ephemeris by using all data and varying 
the period and epoch of periastron. The eccentricity 
and angle of periastron were held fixed for MT91 696, 
CPR2002 A36, and Schulte 3, but varied for MT9f 372. 
Finally, we obtained our best-fit solution by varying the 
remaining variables that depend on both spectroscopic 
and photometric data (i.e., inclination, separation, and 
relative luminosities). Owing to uncertainties on the pri¬ 
mary component temperatures, a small grid of solutions 
was computed for a range of plausible Ti values. When 
a final orbital solution was obtained, we compared the 
results with those obtained using Nightfall, an eclips¬ 
ing binary package created by Rainer Wichmann^^ and 
based on the Wilson & Devinney (1971) code assum¬ 
ing typical Roche geometry. All results from Nightfall 
were consistent with PHOEBE within our final uncertain¬ 
ties. We obtained formal parameter uncertainties via the 
instructions in the PHOEBE manual, where fit parameter 
uncertainties come from the covariance matrices and de¬ 
rived parameter uncertainties are calculated via propa¬ 
gation of error. We note that the solutions of Schulte 3 
and CPR2002 A36 suffer degeneracies stemming from 
probable hot spots and being limited to 1-2 bandpasses. 
Further discussion of these degeneracies are provided in 
Sections 4.5 and 4.6. 

PHOEBE solutions provide computed theoretical abso¬ 
lute bolometric magnitudes. However, because the bolo- 
metric magnitudes represent the components only, rather 
than the components and additional spots, we estimated 
the apparent magnitude for Schulte 3 and CPR2002 A36 
from the model by turning off the spot contribution once 
a solution was finalized. We then used bolometric correc- 

® The O’Connell effect is defined as the unequal out-of-eclipse 
maxima (O’Connell, D. J. K. 1951). Many likely causes have been 
proposed for the phenomenon, but for W Uma systems the proba¬ 
ble cause is a hot spot resulting from a mass stream between the 
two components. 

The Struve-Sahade is defined as the variation in a compo¬ 
nent’s spectral profiles as a function of phase. The equivalent width 
of the lines may or may not vary, and the effect can affect all lines 
or a few. Helium absorption lines are often affected, and the effect 
is generally associated with the secondary. Linder et al. (2007) and 
Bagnuolo et al. (1999) provide a good summary with examples. 

http://www.hs.uni-hamburg.de/DE/Ins/Per 
/Wichmann/Nightfall.html 


tions interpolated from the tables of Martins et al. (2005) 
and the PARSEC isochrones (Bressan et al. 2012) to de¬ 
termine the absolute visual magnitudes. The distance 
to each system may then be derived using the resulting 
distance modulus in conjunction with an estimate of the 
extinction. Ay (§4.2). 

4.2. Extinction Calculations 

We estimate the extinction toward each system from 
the UBVJHK broadband photometry, compiled in Ta¬ 
ble 1 from published sources and from our own UBV pho¬ 
tometry, in conjunction with an adopted spectral energy 
distribution from Kurucz (1992) atmospheres of the ap¬ 
propriate Tf.fj and log(g) for its spectral type. The JHK 
photometry are drawn from the Two-Micron All Sky Sur¬ 
vey (2MASS) Point Source Catalog^®, while UBV mea¬ 
surements come from various independent sources. We 
redden the model atmosphere by either a Cardelli et al. 
(1989, CCM89) or a Fitzpatrick & Massa (2007, FM07) 
prescription for interstellar extinction. Cardelli et al. 
(1989) provide an analytical description of the average 
interstellar extinction law over the ultraviolet, optical 
and near-infrared range, parameterized in terms of two 
variables. Ay and Ry. Fitzpatrick & Massa (2007), by 
contrast, present a table of 328 discrete reddening curves 
toward individual Galactic sources; they highlight the 
uniqueness of each sightline, stressing how generalized 
curves such as Cardelli et al. (1989) are gross averages 
over diverse interstellar dust components. We fit red¬ 
dened photospheres to the data for each star using each 
of these prescriptions. Table 1 the best-fit Ay and Ry 
values. Our approach is more fully described in Vargas- 
Alvarez et al. (2013), who find good agreement between 
these two descriptions of extinction toward the heavily 
reddened Galactic cluster Westerlund 2. They note how 
the non-standard value of i?v=3.8 along that sightline re¬ 
sults in a significantly smaller distance than if the canon¬ 
ical i?v=3.1 were adopted. 

Because these targets are variable at the level of 0.5, 
0.6, 0.3, and 0.5 mag level for MT91 372, MT91 696, 
Schulte 3, and CPR2002 A36 respectively, a robust red¬ 
dening solution ostensibly depends on acquiring UBV 
and JHK datasets at the same orbital phase. The 
2MASS JHK observations were made simultaneously at 
orbital phase, (/)=0.74, 0.84, 0.69, 0.68 for MT91 372, 
MT91 696, CPR2002 A36, and Schulte 3, respectively, 
assuming that the orbital periods have remained con¬ 
stant since the 1998 June 2MASS observations. The 
times of observation are not available in published 
sources for the majority of the UBV photometry, except 
for our new data presented here which we acquired at or¬ 
bital phase (/)=0.32 and 0.47 for MT91 696 and Schulte 3, 
respectively, near JD=2456271.58 early on 2012 Decem¬ 
ber 10 (UT). Assuming that the orbital periods have re¬ 
mained constant and that the broadband colors remain 
unchanged throughout the eclipse cycle (a plausible as¬ 
sumption given the similarity in temperatures between 

Because of the low angular resolution of 2MASS, there may 
be flux contributions from unresolved components. However, only 
in the case of Schulte 3 is there a tertiary bright enough and close 
enough to be considered significant (12% of the system light), and 
even this contribution is inconsequential to the reddening results 
because of the very large extinctions. 
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components in all cases), we correct our UBV photom¬ 
etry to the orbital phase of the 2MASS photometry be¬ 
fore fitting for the optimal reddening solution. These 
corrected values are listed and footnoted as such in Ta¬ 
ble 1. In practice, we find that the reddening toward 
all targets is sufficiently large that adjusting photometry 
to a common orbital phase is a small correction com¬ 
pared to the Ay >5 mag interstellar contributions to the 
broadband magnitudes. Hence, these corrections make 
negligible difference on the final reddening parameters. 

The UBV photometry reported in the literature for a 
given system shows dispersion beyond the formal uncer¬ 
tainties (0.01 - 0.03 mag based on published sources), 
consistent with data being obtained at different orbital 
phases and/or real intrinsic variations in the luminosity 
of the sources. Furthermore, MT91 note the systematic 
differences between their UBV photometry and that of 
Reddish et al. (1966). For purposes of chi-squared min¬ 
imization, we adopt somewhat generous uncertainties of 
0.05, 0.05, 0.03, 0.01, 0.01, 0.01 mag on the UBVJHK 
photometry, respectively. 

Table 1 lists the photometry and best-fit reddening pa¬ 
rameters Ay and i?y for each combination of UBV mea¬ 
surements, paired with the 2MASS JHK photometry, 
for both the CCM89 and FM07 reddening curves. Inspec¬ 
tion of the individual values reveals generally good agree¬ 
ment between the reddening parameters computed using 
different datasets and the same reddening law. For ex¬ 
ample, the extinction results using CCM89 for MT91 372 
are Ay=7.15 and 7.03 mag for the two available UBV 
datasets, giving a mean of 7.09. The results using the 
FM07 reddening curves are Ay=7.01 and 7.13 mag. The 
means of Ay=7.09 and 7.07 mag, respectively, for the 
CCM89 and the FM07 curves are nearly identical. The 
results using CCM89 for MT91 696 are Av=5.82, 6.30, 
and 5.83 with a mean of 5.98. The results using the FM07 
reddening curves are Ay=5.88, 5.94, and 5.75 with a sim¬ 
ilar mean 5.85. CPR2002 A36 and Schulte 3 have mean 
reddenings of Ay=6.74 and 6.07 mag, respectively, using 
the CCM89 reddening curve and Ay=6.46 and 6.02 mag, 
respectively, using the FM07 curves. We conclude that 
there is good agreement between the extinctions derived 
using either the MT91 UBV photometry in conjunc¬ 
tion with 2MASS JHK and our own photometry with 
2MASS. There is also good agreement between the ex¬ 
tinctions derived for MT91 372, MT91 696 and Schulte 3 
(dispersions ^0.1 mag) using the two reddening curves, 
but the dispersion is somewhat larger for CPR2002 A36, 
leading to a larger uncertainty on Ay. 

Finally, we note that the best-fit FM07 reddening 
curves have reduced chi-squared values in the range 1 - 
10 while those for the CCM89 parameterization have 
much larger reduced chi-squared values of 7 - 50. We 
adopt 4 degrees of freedom for the 6 data points in each 
case, although this is not strictly correct for FM07 be¬ 
cause these curves are fitted empirical characterizations 
rather than parameterizations. Nevertheless, the FM07 
series of curves always provides a better fit to the data 
than the average curve of CCM89.^^ The fitting pro¬ 
cess consistently selects curves from the FM07 library 

The Cygnus region may have an anomalous reddening curve 
(Turner 1989) so the CCM89 average parameterization could be 
less applicable along this sightline. 


corresponding to early type stars from NGC 6530, sug¬ 
gesting that this sightline produces reddening similar 
to that of the dust toward Cyg OB2. The mean ex¬ 
tinction and reddening of the FM07 hts for MT91 372 
are 7.07 ± 0.12 mag and 3.37 ± 0.44, respectively. For 
MT91 696, the FM07 fits yield means of Ay = 5.85 ± 
0.09 mag and Ry = 3.30 ± 0.20, and this extinction 
estimate compares favorably to the Ay = 5.84 mag 
found by Massey & Thompson (1991), despite their as¬ 
sumption of Ry = 3.0. The mean of the FM07 fits 
for CPR2002 A36 provide Ay = 6.46 ± 0.23 mag and 
Ry = 3.42 ± 0.19. And finally, for Schulte 3 we find 
Ay = 6.02 ± 0.09 mag and Ry = 3.19 ± 0.23, which is 
slightly smaller than the Ay = 6.24 mag found by Torres- 
Dodgen et al. (1991). Torres-Dodgen et al. (1991) deter¬ 
mined i?v=3.04±0.09 for dozens of stars in the Cyg OB2 
vicinity. Our average value using the CCM89 redden¬ 
ing curve is i?y~3.05±0.05, whereas our average value 
using FM07 is i?v~3.3±0.1, somewhat higher but con¬ 
sistent with values found throughout the region. Here¬ 
after, we adopt the best-fit FM07 extinction fits. We 
note that adoption of the CCM89 reddening law would 
systematically lower the derived distance modulus for 
CPR2002 A36 by about 0.28 mag, but it will have a very 
small impact on the distances for MT91 372, MT91 696 
and Schulte 3. 

4.3. MTH 372 

Of the four eclipsing systems analyzed, MT91 372 is 
one of the most ideal for determining distance. Figure 4 
shows the velocity curve data and best-ht model (top), 
light curve and best-fit model (middle), and light curve 
residuals (bottom). The components are detached with 
clear, deep eclipses having Am ~ 0.2 for the secondary 
eclipse and Am ~ 0.5 for the primary eclipse in V. Fig¬ 
ure 5 shows disentangled spectra for each component. 
The primary star displays strong He I A5876 lines rela¬ 
tive to Ha, indicative of a hot star in the range BOV - 
B2V. While He I is present in the secondary, this ratio is 
smaller, indicating a later spectral type near mid-B. 

Because this system is among the faintest and red¬ 
dest objects in the Cygnus OB2 Radial Velocity Survey 
{y = 14.96 at its brightest and Ay—7 mag), the SNR 
of the Survey’s WIYN spectra are insufficient to use the 
well-calibrated FW 448 i/£’W 447 i ratios to measure com¬ 
ponent temperatures. We determine a 2cr upper limit of 
Ti < 28, 000 K for the primary based on an upper limit 
for the ratio EW^^^ll/< 0.10 (calibrated using 
Figure 3) in a combination of all available WIRO spectra 
taken near the same quadrature phase.The available 
spectra provide no strong lower limits on the effective 
temperature of the primary. The presence of He I A5876 
with an EW of 0.60±0.02 A in the separated primary 
star spectrum in Figure 5 confirms that MT91 372a is 
a hot star with Ti > 18, 000 K based on the calibration 

The measured (indicated by subscript m) equivalent widths 
of lines in a composite spectrum of two stars, EW^^i and EWTn, 2 -, 
yield the true (indicated by subscript t) equivalent widths via 
the ratio of continuum light at that wavelength, I = C\IC 2 ^ as 
EWyi = EWmA^ + 0/^ and EWt ,2 = EW^M^ + 0 (Hilditch 
2001, eqn. 5.98). However, assuming that the continuum light ra¬ 
tios of the two very hot stars at 5876A is the same as at 541lA, a 
direct measurement of the ratio EWs^ii/EW^sTe for each star is 
essentially independent of 1. 
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TABLE 1 

Ay AND Ity ESTIMATES 


Quantity 

MT91 372 

MT91 696 

CPR2002 A36 

Schulte 3 

u 

18.01“* 

14.44“* 

13.73“* 

12.31“* 

B 

17.14“* 

13.97“* 

12.78“* 

11.85“* 

V 

14.97“* 

12.32“* 

11.36“* 

10.09“* 

J 

12.51** 

8.53** 

7.19** 

6.50** 

H 

9.92** 

8.14** 

6.65** 

6.00** 

K 

9.60** 

7.89** 

6.36** 

5.74** 

Av(CCM89) 

7.15 

5.82 

6.63 

6.16 

i?v(CCM89) 

3.03 

2.90 

3.11 

3.11 

Av(FM07) 

7.01 

5.88 

6.56 

5.92 

itvlFMO?) 

2.93 

3.06 

3.31 

2.93 

U 

17.19 

14.20“* 

13.73“* 

12.30“= 

B 

16.52 

13.87“* 

12.81“* 

11.83® 

V 

15.04 

12.36“* 

11.27“* 

10.22® 

Av(CCM89) 

7.03 

6.30 

6.96 

5.93 

ii:v(CCM89) 

2.95 

3.11 

2.92 

2.92 

Av(FM07) 

7.13 

5.94 

6.46 

6.07 

A;v(FM07) 

3.81 

3.42 

3.31 

3.30 

U 


14.00* 

13.19* 

12.12* 

B 


13.69* 

12.77* 

11.77* 

V 


12.20* 

11.15* 

10.22* 

Av(CCM89) 


5.83 

6.63 

6.13 

ii:v(CCM89) 


2.98 

3.02 

3.02 

Av(FM07) 


5.75 

6.36 

6.09 

ii:v(FM07) 


3.42 

3.64 

3.36 

Avg. Ay (CCM89) 

7.09±0.12 

5.98±0.27 

6.74±0.19 

6.07±0.12 

Avg. By (CCM89) 

2.99±0.08 

3.00±0.10 

3.09±0.10 

3.02±0.09 

Avg. Ay (FM07) 

7.07±0.12 

5.85±0.09 

6.46±0.23 

6.02±0.09 

Avg. Ry (FM07) 

3.37±0.44 

3.30±0.20 

3.42±0.19 

3.19±0.23 


^ Massey & Thompson (1991) 

^ 2MASS; data acquired at orbital phase (f>—0.74, 0.84, 0.69, 0.68 for MT91 372, 
MT91 696, CPR2002 A36, and Schulte 3, respectively, assuming that the orbital 
periods have remained constant since the 1998 2MASS observation. 

^ Reddish et al. (1966) 

Tycho-2 H0g et al. (2000) 

^ Nicolet (1978) 

^ This work; obtained at orbital phase 0—0.32,0.49,0.47 for MT91 696, 
CPR2002 A36, and Schulte 3, respectively We have corrected our measured UBV 
magnitudes to the same orbital phase as the 2MASS JHK measurements using 
the derived orbital solution under the assumption of constant orbital periods and 
constant broadband color. 


between Tf,ff and He I A5876 presented in Figure 3 of 
Kobulnicky et al. (2012). The ratio — 5 

provides similar loose constraints. 

In an attempt to provide a more secure lower limit on 
the effective temperature of the primary, MT91 372 was 
observed at WIRO on 2013 August 14 and August 29 for 
3 hours (12 x 900 s) each using a 2000 1 mm“^ grating 
in first order over the wavelength range 6900 - 8000 A. 
The dates of observation were chosen to cover portions 
of the orbital phase (0.09 and 0.85, respectively) where 
the temperature-sensitive far-red lines of He I A7065, 
O I AA7772,7774,7775 A, and Mg II A7896 fell at lo¬ 
cations minimally affected by atmospheric absorption, 
interstellar bands, and night sky emission lines. No or¬ 
der blocking filter was used as the atmosphere, detector, 
and heavy interstellar extinction {Ay = 7 mag) effec¬ 
tively serve to block wavelengths shorter than 4000 A. 
Data reduction followed standard procedures described 
in Section 2. Particular care was taken with the con¬ 
tinuum normalization of the spectra to ensure that line- 
free portions were used to fit a low-order polynomial to 
the spectrum. Regions of strong atmospheric absorption 
near 7150 - 7400 A and 7600 - 7700 A were excluded from 


the fit. From these spectra we detect He I A7065 at the 
expected level for stars between 20,000 K and 28,000 K, 
but O I AA7772,7774,7775 A, and Mg II A7896 lines are 
not detected to upper limits of ~0.12 A. While this limit 
does not place useful constraints on the temperature on 
the basis of the expected Mg II A7896 equivalent widths 
(0.05 A or less for stars hotter than 20,000 K, based on 
Tlusty models), the limits are useful in the case of the 
O I AA7772,7774,7775 blend which ought to have a com¬ 
bined EW>0.2 A for stars cooler than 20,000 K. 

We conclude that 20,000 K is a secure lower limit on Ti. 
Accordingly, we adopt a temperature of 24,000±3,000 K 
for the primary component, MT91 372a. Given that the 
system is clearly detached with a small separation, the 
components are not likely filling their Roche lobes. Based 
on the absence of emission features we assume that the 
components are un-evolved. The joint light curve and 
velocity curve analysis to follow confirms that the com¬ 
ponents have radii consistent with main-sequence, early 
B stars. 

Radial velocities obtained from WIRO-Longslit spectra 
between 2008 and 2013 (12 epochs for the primary, 11 
epochs for the secondary) produce power spectra with 
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Fig. 4.— Radial velocity curve (top), light curve (middle), and 
light curve residuals (bottom) for MT91 372. Solid and unfilled 
points (top) denote the primary and secondary radial velocities, 
respectively. Solid curves in each panel show the best-fit PHOEBE 
solutions given the adopted parameters, as described in the text. 
The steps at phases of 0.5, 1.0, and 1.5 are a result of the Rossiter- 
McLaughlin effect (Rossiter 1924; McLaughlin 1924) 



5400 5600 5800 6000 6200 6400 6600 

Wavelength (A) 


Fig. 5. — Normalized component spectra for MT91 372a and 
MT91 372b using the Gonzalez & Levato (2006) technique for dis¬ 
entanglement. The spectra are shown with arbitrary vertical off¬ 
sets. 


two strong signals, the likely period at ~2.218 d and the 
\—v alias at ^1.821 d. The BSCSP solution also converged 
on a 2.218 d period. 

Our photometric and spectroscopic data, adopted Ti, 
and starting orbital parameter guesses were used to com¬ 
pute the best-fit PHOEBE orbital solution displayed in Fig¬ 
ure 4. The best period for the combined datasets is P = 
2.227490 ± 0.000005 days, and the light curve requires a 
non-zero, but low, eccentricity of e = 0.044 ± 0.002. The 
best PHOEBE solution provides a slightly lower mass ratio 
oiq = 0.44 ± 0.01 (BSCSP yields q = 0.50 ± 0.02). Ta¬ 
ble 2 lists the best-fit orbital parameters, including the 
orbital period (P), the date of periastron (Po)j eccen¬ 
tricity (e), orbital inclination (i), semi-major axis (a), 
systemic velocity ( 7 ), and mass ratio {q). Table 2 also 
lists physical parameters for the system derived from the 
modeling process, including the secondary temperature 
(P 2 )) the radii of both components (Pi and P 2 ), the 
masses (Mi and M 2 ), and the bolometric magnitudes 
{MboIi and Mson)- Not surprisingly, given the depth 
of the eclipses, MT91 372 is nearly edge-on with an in¬ 
clination of i = 86 ± 1°. This results in computed com¬ 
ponent masses of 13.0 ± 0.7 Mq and 5.7 ± 0.4 Mq and 
radii of 4.94 ± 0.04 Rq and 3.34 ± 0.06 R©. Coupled 
with the estimated primary temperature and the result¬ 
ing best fit secondary temperature (P 2 = 16,600±2,300), 
the best-fit spectral types are BIV and B5V (Drilling & 
Landolt 2000), agreeing well with our initial estimates. 
The hnal rows of Table 2 list the adopted extinction 
{Ay=7.07±0.12), and the observed out-of-eclipse appar¬ 
ent magnitude of the system {my = 14.96), the distance 
modulus (D.M.), and the resulting distance (d) in kpc. 
Given the derived bolometric magnitudes (—4.87 mag 
and —2.41 mag for the primary and secondary, respec¬ 
tively), the distance modulus may be computed in the 
standard manner, 

D.M. = mv — My — Ay 

= my + 2.h \og{lO-^^Bon-BC^)/2.5 

-|- 2Q-(*fBoi2-BC'2)/2.5^ _ 

The resulting distance modulus of 10.62 ±0.28 equates 
to a distance of 1.33 ± 0.17 kpc. We note here that the 
uncertainties on MboIi and Mbo 12 are coupled, not in¬ 
dependent, since T 2 is derived from Ti via the ratio of 
eclipse depths. Furthermore, there is an anti-correlation 
between bolometric magnitudes (i.e., temperatures) and 
the bolometric corrections (B.C.), meaning that an er¬ 
ror in MboI is partially canceled by a change in B.C. in 
the opposite direction. Consequently, the derived My 
used to determine the distance modulus is less sensitive 
than one might ordinarily assume to the adopted primary 
effective temperature, Ti. Nevertheless, the relatively 
large uncertainty on the temperature of MT91 372A is 
the dominant source of uncertainty on the derived dis¬ 
tance for this system. 

4.4. MT91 696 

When corrected for third light contributions, MT91 
696 is also a good object for determining the distance 
to Cyg OB2. As a relatively bright {my = 12.05 mag), 
short-period system (P = 1.46918 d), it contains com¬ 
ponents close in spectral type (09.5V and B0.5V) and 
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luminosity (Kiminki et al. 2012), the system is nearly 
edge-on with eclipse depths of ^ 0.64 mag. Figure 6 dis¬ 
plays the radial velocity curves (top), light curve (mid¬ 
dle) , and light curve residuals (bottom). The solid curves 
represent the best-fit PHOEBE model discussed below. The 
smooth, continuous variation in the light curve indicates 
that the components are near contact. Velocity and pho¬ 
tometric uncertainties are small compared to the varia¬ 
tions. 



Fig. 6.— Radial velocity curve (top), light curve (middle), and 
light curve residuals (bottom) for MT91 696. Solid and unHlled 
points (top) denote the primary and secondary radial velocities, 
respectively. Solid curves in each panel show the best-Ht PHOEBE 
solutions given the adopted parameters, as described in the text. 

To estimate the component temperatures, we utilized 
21 epochs of WIRO-Longslit spectra taken over nine 
nights between 2008 August 20 and 2012 October 21 to 
disentangle the component spectra (shown in Figure 7) 
and obtain the cross-correlation radial velocities listed 
in Table A.l. In the disentangled primary spectrum, the 
EW ratio of temperature-sensitive lines EW^/^wj 
is 0.45 ± 0.05. We estimate the ratio for the secondary 
to be FTF 54 ii/ifIF 5876 = 0.25 ± 0.05. We also mea¬ 
sured the equivalent width ratio using a composite of 
four spectra obtained near a single quadrature phase 
((/> = 0.75). The measured ratios were higher using this 
method, 0.76 ± 0.07 and 0.53 ± 0.07 for the primary and 
secondary respectively. We adopted the mean of these 
ratios, 0.61 ± 0.09 for the primary and 0.44 ± 0.09 for 
secondary. Using the relationship shown in Figure 3, this 
yields Ti = 32,000 ± 600 K and T 2 = 29,800 ± 700 K. 
The primary temperature is similar to the 31,884 K im- 



Fig. 7. — Normalized component spectra for MT91 696A and 
MT91 696B using the Gonzalez & Levato (2006) technique for dis¬ 
entanglement. The secondary spectrum is shown with an arbitrary 
vertical offset for clarity. 


plied by the observational scale of Martins et al. (2005) 
for an 09.5V. Accordingly, we adopt T\ = 32,000 ± 600 
as a hxed parameter in the PHOEBE modeling but leave 
Ti as a free parameter to be fit. 

Before considering the third light contribution from 
MT91 696C, we computed the best-fit PHOEBE model to 
the photometry that includes this extra component. We 
then tested various orbital configurations, and each indi¬ 
cated that the system is just barely detached. Once the 
code converged on the best-fit solution, we added increas¬ 
ingly larger contributions from third light until PHOEBE 
could no longer converge on a solution and the primary 
eclipse depth exceeded 0.75 magnitudes. The highest 
contribution allowed was 14% of maximum out-of-eclipse 
light (i.e., 14% of the V=12.051 mag system). This in¬ 
dicates that the tertiary cannot be brighter than about 
V=14.18 mag. Our ephemeris indicates that the EIST 
FGS measurement was taken at a phase near maximum 
light. The 0.94±0.40 mag offset between the (A-PB) bi¬ 
nary and the brightest astrometric companion, as re¬ 
ported by Caballero-Nieves et al. (2014), results in a 
nominal contribution of 22^^° % — consistent with the 
maximum of 14% allowed by the light curve. We there¬ 
fore adopt a contribution of 14% for component c. 

Table 2 displays the resulting system parameters for 
the best solution to the MT91 696 system after includ¬ 
ing a constant 14% third light contribution within the 
PHOEBE model and fixed zero eccentricity (when allowing 
the eccentricity to vary, BSCSP and PHOEBE both con¬ 
verged on the same zero eccentricity solution). The final 
orbital period obtained with all radial velocity and light 
curve information is in excellent agreement with the pe¬ 
riod reported by Souza et al. (2014), indicating there has 
been negligible change in the period over six years. Addi¬ 
tionally, while the computed radii (i?i = 6.36 ± 0.03 Rq 
and i ?2 = 5.72 ± 0.02 Rq) are slightly smaller than the 
theoretical values of 7.3±0.1 R© (09.5V) and 6.6±0.1 R© 
(B0.5V), the computed masses (Mi = 15.0±0.2 Mq and 
M 2 = 12.1±0.8 M©) are in good agreement with the the- 
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TABLE 2 

Joint Light and Velocity Curve Analysis for MT91 372, MT91 696, CPR2002 A36, and 

Schulte 3 



MT91 372 

MT91 696 

CPR2002 A36 

Schulte 3 

P (days) 

2.227582±0.000005 

1.469179±0.000001 

4.6749±0.0006 

4.74591±0.00005 

To (HJD-2,400,000) 

56100.854±0.001 

54731.6401±0.0005 

55725.8±0.2 

53996.093±0.006 

e 

0.044±0.002 

0 (fixed) 

0 (fixed) 

0 (fixed) 

A (Rq) 

19.0±0.3 

16.3±0.2 

38.7±0.8 

36.0±0.9 

q 

0.44±0.01 

0.80±0.02 

0.62±0.02 

0.42±0.03 

Vo (km s“^) 

-7±2 

1.1±2.7 

-29±3 

-47±3 

i (degrees) 

86±1 

86.5±1.5 

73±1 

59.2±0.3 

3rd light 


14% (fixed) 


12% (fixed) 

log(9) (P) 

4.2 

4.0 

3.4 

3.3 

log(9) (S) 

4.1 

4.0 

3.4 

3.3 

Ti (K) 

24,000±3,000 

32,000±600 

30,000±2,000 

38,500 

T 2 (K) 

16,600±2,300 

30,940±1,050 

20,600±600 

31,300±700 

Ri (R©) 

4.94±0.04 

6.36±0.03 

16.4±0.1 

16.5±0.1 

R 2 (R©) 

3.34±0.06 

5.72±0.02 

13.1±0.1 

ll.OitO.l 

Ml (M©) 

13.0±0.7 

15.0±0.2 

22±2 

20±2 

M 2 (M©) 

5.7±0.4 

12.1±0.8 

14±1 

8.2±0.9 

Spectral Type 1 

BIV 

09.5V 

09 - 09.5IIL 

06V 

Spectral Type 2 

B5V 

B0.5V 

09.5 - BOIV: 

09III - 0911 

MboH (mag) 

-4.87±0.60 

-6.66±0.14 

-8.44±0.30 

-9.54±0.54 

Mbo 12 (mag) 

-2.41±0.68 

-6.29±0.20 

-6.32±0.30 

-7.76±0.39 

B.C.i (mag) 

-2.38±0.33 

-3.05±0.09 

-2.88±0.16 

-3.59±0.37 

B.C .2 (mag) 

-1.45±0.38 

-2.96±0.10 

-2.01±0.17 

-2.99±0.15 

Mv,i (mag) 

-2.49±0.27 

-3.61±0.04 

-5.56±0.15 

-5.95±0.17 

Mv ,2 (mag) 

-0.96±0.30 

-3.33±0.10 

-4.32±0.13 

-4.77±0.25 

Mv (mag) 

-2.73±0.28 

-4.24±0.07 

-5.86±0.14 

-6.26±0.19 

Av (mag) 

7.07±0.12 

5.85±0.09 

6.46±0.23 

6.02±0.09 

mv (mag) 

14.96±0.02 

12.22±0.01=" 

11.40±0.0l'" 

10.36±0.03‘‘’'" 

RMSv (mag) 

0.023 

0.007 

0.015 

0.013 

RMS NSVS (mag) 



0.035 

0.043 

RMSrvi (kms-l) 

10.2 

12.9 

25.2 

32.9 

RMSrv 2 (km s-l) 

14.3 

19.3 

20.5 

37.6 

D.M. (mag) 

10.62±0.28 

10.61±0.10 

10.80±0.27 

10.60±0.21 

Dist. (kpc) 

1.33±0.17 

1.32±0.07 

1.44±0.18 

1.32±0.13 


^ component A+B only 
^ spot contribution removed 


oretical values of 16.0 ±0.5 Mq and 12.9 ±0.4 Mq (Mar¬ 
tins et al. 2005). Finally, the resulting secondary tem¬ 
perature, ±2 = 30,940 ± 1050 K also agrees well with the 
direct T 2 measurement from the He II A5411/He I A5876 
line ratio determined above. 

The light curve residuals in the lower panel of Figure 6 
show small <1% systematic variations between phase 0 
and 0.25, suggesting minor deficiencies in the adopted 
model, but additional high-quality photometry will be 
needed to further constrain the model parameters. The 
computed distance modulus is 10.61±0.10 mag and cor¬ 
responds to a distance of d=1.32±0.07 kpc. For this 
system, the largest source of uncertainty on the distance 
is the uncertainty on the correction for interstellar ex¬ 
tinction. 

4.5. CPR2002 A36 

Ideally, eclipsing binary distances would make use of 
detached systems which have simple geometries, i.e., the 
stars are spherical and the light curves are unaffected by 
complications involving mass loss, hot spots produced 
by radiation or mass transfer streams, and variations in 
the surface potentials of stars approaching or undergo¬ 
ing Roche lobe overflow. Unfortunately, four of the six 
known eclipsing SB2 systems in the direction of Cyg OB2 
contain evolved components and are probable contact or 
over-contact systems. Therefore, we proceed to model 


the CPR2002 A36 and Schulte 3 eclipsing SB2 systems 
while cautioning that the interpretation of the velocity 
curves and light curves is subject to larger uncertainties. 

Figure 8 shows the best-fit light curve data and model 
(top), velocity curve data and model (middle), and light 
curve residuals (bottom) for CPR2002 A36. The light 
curve is characteristic of an over-contact system with a 
difference in eclipse depths that indicates a significant 
difference in component effective temperatures. Addi¬ 
tionally, we see an O’Connell effect that results in an 
^0.04 mag difference in maxima. There is also con¬ 
siderable dispersion in the NSVS magnitude at nearly 
all phases, exceeding the nominal 0.01 mag uncertain¬ 
ties and suggestive of a high level of intrinsic temporal 
variability. Finally, because of the over-contact nature 
of the system, CPR2002 A36 shows strong, varying Ha 
emission and a Struve-Sahade effect in the He I lines. 

The emission and Struve-Sahade effect are mostly re¬ 
sponsible for the disagreement between the observed and 
calculated radial velocities. Emission partially fills the 
spectral lines and skews the He I line centers and is most 
signihcant for phases ~0.8 - 1.1 when the secondary is 
heading into eclipse. This emission is easiest to see at Ha, 
shown in the panel second from right in Figure 9 (colored 
white in this inverse greyscale image), a phased spectral 
sequence using all of the WIRO-Longslit spectra folded 
with the final solution period. The figure also shows 
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Fig. 8. — Radial velocity curve (top), light curve (middle), and 
light curve residuals (bottom) for CPR2002 A36. Solid and unfilled 
points (top) denote the primary and secondary radial velocities, 
respectively. Solid curves in each panel show the best-fit PHOEBE 
solutions given the adopted parameters, as described in the text. 


four additional spectral lines in the AA 5400 - 5700 A 
range: He II A5411, O III A5592, He I A5876, and Ha, 
and He I A6678. The orange and green curves represent 
the final primary and secondary solutions, respectively. 
From a phase of ^0.25 - 0.6, a singular emission fea¬ 
ture is evident and appears to be associated with the 
redshifted secondary. At phase 0.85 - 1.1, there appear 
to be multiple emission components with the dominant 
one still associated with the secondary. At He I A5876 
(middle panel), the emission results in the best Gaussian 
fit being centered at a longer wavelength than the true 
center. Because the primary component’s He I absorp¬ 
tion is weaker and shifted less, the He I profiles are more 
distorted at these phases, and therefore, computed radial 
velocities are more skewed. The emission, Struve-Sahade 
effect, and O’Connell effect may be adequately explained 
by a hot spot or a stream of optically thin ionized gas 
near or trailing the secondary. 

An advantage to creating a phased spectral image 
is that weak lines that would normally be difficult to 
identify amongst the continuum noise and various in- 
terstellar/DIB lines are more easily identified and as¬ 
sociated with the correct component. For instance, 
O III A5592 is difficult to identify with certainty in 
many of CPR2002 A36’s spectra, but its presence is 
easily revealed in both components in Figure 9. The 
full phased spectral image also reveals moderately weak 
C IV AA580I, 5812 absorption and very weak C III A5696 



Fig. 9.— A sequence of spectra for CPR2002 A36 ordered by 
orbital phase, showing five spectral lines with the spectroscopic 
orbital solution overplotted in orange for the primary and green 
for the secondary. 


emission in the primary, the presence and strength of 
which, with He II A5411, indicate the primary may be 
slightly hotter than the BO estimated in Kiminki et al. 
(2009), such as 09 - 09.5 (Walborn 1980). Given the 
strength of the metal lines in these spectra and of the 
Si IV and C III lines shown in Figure 4 of Kiminki et al. 
(2009), we revise our estimate of the primary’s luminos¬ 
ity class to HI and adopt a full spectral classification of 
09 - 09.5HI. For the secondary, the full phased spectral 
image revealed the presence of weak N II AA5667, 5670, 
5680. Given these, the weak He II A5411, and the very 
weak O III A5592, we revise and adopt the spectral type 
of 09.5 - BOIV: for the secondary. 

We estimate the effective temperatures for this system 
in a similar manner to MT91 696 via the ratio of he¬ 
lium equivalent widths. We measured equivalent widths 
at both quadrature phases, obtaining AIF 5411 /AW 5876 = 
0.52±0.09 and EW^/^wj = 0.48±0.04 for the pri¬ 
mary and secondary, respectively. As both components 
of CPR2002 A36 have been classified as likely evolved, we 
use the dashed line in Figure 3 to estimate the tempera¬ 
tures as Ti=29,400 K and T2=29,200 K with uncertain¬ 
ties of approximately 500 K. These estimates are in rough 
agreement with the theoretical temperatures for an 09HI 
(31,000±800 K) and 09.5IV (30,800±500), but disagree 
with the temperature difference indicated by the differ¬ 
ence in eclipse depths. However, considering the lower 
signal-to-noise of He II A541I in many of the spectra and 
the fact that the Struve-Sahade effect raises the inherent 
uncertainty in equivalent width measurements of He I, 
and given the theoretical temperatures (Martins et al. 
2005) for the estimated spectral types above, we adopt 
30, 000 ± 2, 000 K for Ti and estimate 29, 000 ± 2, 000 K 
for T 2 . 

Our nominal PHOEBE model is a double-contact con¬ 
figuration with a hot spot on the trailing side of the 
secondary (130° from the point of contact on the equa- 
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tor). In this scenario, the visibility of the spot coin¬ 
cides well with the presence of emission in Figure 9. The 
spot may indicate the transfer of mass from the primary. 
The final solution is fit with a spectroscopic period of 
P = 4.6749 ± 0.0006 d and eccentricity fixed at zero. We 
were unable to reconcile a phase shift of 0.157 (17.6 hr) 
between the NSVS photometry and our radial velocity 
data with PHOEBE. Therefore, we adopted the ephemeris 
obtained from BSCSP and applied an arbitrary 0.1565 
phase shift to the NSVS photometric data. Similarly, 
we applied an arbitrary 0.0506 phase shift to the RBO 
V-band photometry as well. It is possible that the sys¬ 
tem has a changing period that we have not accounted 
for, but its impact is negligible on the results presented 
here. Additionally, we find that the V-band photometry 
and the radial velocity curves are better fit with a slightly 
higher inclination than the NSVS photometry {i = 72.9° 
versus i = 69.6°). Because NSVS photometry is filter¬ 
less and only approximate to Johnson-R band (Wozniak 
et al. 2004), we adopt the solution with the higher in¬ 
clination. This also has negligible impact on the results 
presented here. In Figure 8, we show the NSVS pho¬ 
tometry instead of the V-band photometry because its 
full phase coverage better illustrates the solution’s over¬ 
all agreement with the observed light curves. 

The final solution yields a secondary temperature of 
T2=20,600±600 K, contrasting significantly with the 
29, 000 ± 2,000 K estimated by helium equivalent width 
ratios, but relative to Ti, it is in better agreement with 
the temperature difference indicated from the eclipse 
depths. No combination of spot parameters and overall 
model allowed for a higher secondary temperature, but 
we admit that fitting the spot manually limits the pa¬ 
rameter space explored. The computed masses Mi=22± 
2 Mq and M 2 =14 ± 1 Mq are in good agreement with 
the expected values from Martins et al. (2005). However, 
the radii, i?i=16.4 ± 0.1 R© and i?2=13.1 ± 0.1 R©, are 
larger than the expected values and more characteristic 
of 0911 and 09.5III. However, given the uncertainties on 
the spectral types, and the fact that each component is 
filling its Roche lobe, these values are not entirely surpris¬ 
ing. Using the computed bolometric magnitudes in Ta¬ 
ble 2, we compute a distance for this system of 1.44±0.18 
kpc, where the uncertainty on distance is dominated by 
the relatively large uncertainty on interstellar extinction 
((tai/= 0.23 mag) and the uncertainty on primary tem¬ 
perature (ctt =2000 K). Had we adopted the CCM89 
reddening law for CPR A36 instead of FM07, the result¬ 
ing distance would drop to 1.27 kpc. 

4.6. Schulte 3 

Schulte 3 shows many similarities to CPR2002 A36, in¬ 
cluding highly variable Ha emission, a Struve-Sahade ef¬ 
fect, an O’Connell effect, and a similar period (4.74591 ± 
0.00005 d). Figure 10 shows the velocity curve data and 
model (top), light curve data and model (middle), and 
the light curve residuals (bottom). In the top panel, 
squares represent WIRO data, triangles represent WIYN 
data, and stars represent WIRO data excluded when de¬ 
termining the BSCSP and PHOEBE solutions. We excluded 
these data because of the distortion of the He I A5876 
profiles near quadrature. Like CPR2002 A36, this dis¬ 
tortion, the Struve-Sahade effect, more strongly affects 
the primary radial velocities and is more easily seen in 


Ha, shown in the bottom panels of Figure 11 (colored 
white in this inverse greyscale image). The figure is 
a phased spectral sequence using all of the Hydra and 
WIRO-Longslit spectra folded with the final solution pe¬ 
riod, where the orange and green lines represent the pri¬ 
mary and secondary spectroscopic solutions plotted over 
various spectral lines, respectively. The emission in Ha 
and He I appears at slightly different phases than with 
CPR2002-A36, indicating a different spot geometry and 
location. In the Figure, Ha shows a singular emission fea¬ 
ture, which appears to be associated with the redshifted 
secondary, at phases ~0.25 - 0.5. At phases 0.6 - 1.0, it 
shows multiple emission components with the dominant 
one still associated with the secondary. This emission 
is also present in He I A5876, but it is more difficult to 
see in the Figure owing to its weaker presence in the 
wings around this line. Regardless, the primary compo¬ 
nent’s He I radial velocities are skewed at these phases 
for the same reasons given for CPR2002 A36. We ob¬ 
tained 20 epochs of WIYN/Hydra spectra and 24 epochs 
of WIRO-Longslit spectra. Radial velocities were mea¬ 
sured for 12 and 19 of these spectra, respectively. How¬ 
ever, because of the distorted line profiles we utilized only 
WIYN radial velocities for the primary and excluded all 
secondary He I A5876 radial velocities near (f) = 0.25. 
The quality of the final solution is demonstrated with the 
good agreement between the radial velocity curves and 
the spectroscopic solutions in Figure 11. In addition to 
Ha in the bottom-right panel, the top-right panel shows 
Si IV 4089 A, N III 4097 A, RS 4101 A, and Si IV 4116 A, 
and the middle right panel shows He I 5876 A. We pro¬ 
vide all radial velocities in Table A.l, with a note to 
designate the observations excluded from fitting. The 
table also includes the observed minus calculated values 
of the excluded points for reference. 

We also utilized the WIYN spectra, after removing 
contributions from interstellar Ca II and the diffuse inter¬ 
stellar band (DIB) at A 4428 A, to disentangle the com¬ 
ponent spectra. Because the variability of the He I line 
depths more strongly affects the primary than the sec¬ 
ondary, He I is underrepresented in the disentangled pri¬ 
mary spectrum. The spectrum shows EW 4200 /AW 4026 > 
1, indicating a temperature class of early to mid-0. How¬ 
ever, taking into account that He I is underrepresented, 
the ratio may be closer to unity, indicating a tempera¬ 
ture class near 06 - 07. Additionally, the weak N III 
with absent O II, N II, and Si IV, indicates a likely lumi¬ 
nosity class of V. The disentangled secondary spectrum 
shows EW 4200 /£’ 1^^4144 near unity, indicating a temper¬ 
ature class near 09, and the strength of Si IV and N III 
lines relative to the neighboring He I lines indicates an 
evolved luminosity class, likely HI - H. The disentan¬ 
gled component spectra are shown at top and third from 
top in Figure 12. The spectra show good agreement 
with the spectral standards of HD101190 (06V; second 
from top) and HD148546 (09Ia; bottom) from the Wal- 
born & Fitzpatrick (1990) digital atlas. We note that 
HD101190 has a luminosity of la rather than HI, result¬ 
ing in stronger metal lines, but it is the closest standard 
in the library. Interestingly, the secondary luminosity 
class suggests that the present secondary may have been 
the more massive component at one point. 

As with the other systems, we measured the effective 
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Fig. 10. — Radial velocity curve (top), light curve (middle), 
and light curve residuals (bottom) for Schulte 3. Solid and un¬ 
filled points (top) denote the primary and secondary radial ve¬ 
locities, respectively. Squares represent observations with WIRO- 
Longslit, triangles represent observations taken with WIYN-Hydra, 
and stars represent the observations excluded from fitting. Solid 
curves in each panel show the best-fit PHOEBE solutions given the 
adopted parameters, as described in the text. 

temperatures using the ratio of helium line equivalent 
widths in spectra obtained near quadrature. We find 
that EW^iii /is 1.06±0.06 for the primary and 
0.95±0.04 for the secondary. Here, the Struve-Sahade 
effect either does not impact our equivalent width mea¬ 
surements or it impacts them equally, as the variation in 
He I line profile associated with the primary appears to 
be limited to line depth and width changes only; the over¬ 
all equivalent width does not change between quadrature 
phases. With these equivalent width ratios and the lu¬ 
minosity class estimates, the relations provided in Fig¬ 
ure 3 yield temperatures of Ti = 34, 500 ± 500 K and 
T 2 = 32,000 ± 500 K for the primary and secondary, 
respectively. The secondary estimate is nearly in agree¬ 
ment with the theoretical value of 31, 300 ± 600 K from 
Martins et al. (2005). However, the primary estimate is 
lower than the theoretical value for an 06V by approxi¬ 
mately 4,000 K, which makes the temperature more ap¬ 
propriate for a slightly later type of 07 (Martins et al. 
2005). Given that we cannot fully quantify the effect that 
emission plays with the equivalent width ratios, we also 
consider that the effective temperature may be as high 
as 38,500 K for the primary. Therefore, we computed 
a grid of solutions for primary temperatures 34,000 K - 
38,500 K. In these solutions, the secondary temperatures 
from ~28,000 K - 35,500 K depending on spot parame- 


Fig. 11.— A sequence of spectra for Schulte 3 ordered by or¬ 
bital phase, showing nearly the entire wavelength range covered 
by our WIYN and WIRO observations. Right panels illustrate the 
quality of the overall solution with spectroscopic orbital solution 
overplotted in orange for the primary and green for the secondary. 

ters and model used. 

The light curve for Schulte 3 is unambiguously charac¬ 
teristic of a contact/overcontact system. It exhibits rela¬ 
tively small eclipse depths of ~0.23 mag and ~0.I7 mag 
in the raw light curve, suggesting a moderate inclina¬ 
tion. Similarly to CPR2002 A36, the light curve also 
displays an O’Connell effect (~0.02 - 0.03 mag) with 
a brighter primary maximum and exhibits real photo¬ 
metric variability in excess of the nominal uncertainties. 
Photometry obtained at similar orbital phase across dif¬ 
ferent orbits exhibit offsets of up to 0.02 mag, signifying 
intrinsic variability. This is probably a triple system, and 
therefore the light curve (and to some minor extent, the 
velocity curves) is further complicated by the presence of 
Schulte 3C which contributes ~I2% to the V-band light. 
We used PHOEBE to account for the third light contribu¬ 
tion. 

We explored various combinations of temperatures for 
the primary, secondary, and spot, various spot configu¬ 
rations, and both W-Uma and double-contact configura¬ 
tion models for Schulte 3. Because the emission features 
appear to be associated mainly with the secondary, we 
placed the relatively large (49 - 60°) hot spot on that 
component, with a center between 24° and 45° from 
the northern pole and between 20° and 35° counter¬ 
clockwise from the point of contact. The adopted lo¬ 
cation of this hotspot coincides well with the presence of 
emission seen in Ha. In the end, we adopted the solution 
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Fig. 12.— Normalized component spectra for Schulte 3A and 
Schulte 3B using the Gonzalez Sz Levato (2006) technique for dis¬ 
entanglement. Comparison standard spectra from the atlas of Ja¬ 
coby et al. (1984) appear below each component spectrum. The 
spectra are shown with arbitrary vertical offsets. 

with Ti = 38,500 K and T 2 = 31,300 K from Martins 
et al. (2005) because our confidence in the spectral types 
exceeded that of the compromised equivalent width ra¬ 
tios. Additionally, the spot placement for this solution 
is nearer to the point of contact (45° from the northern 
pole and centered 35° in the direction of rotation) than 
the other models, which makes sense if the primary is 
heating the secondary. This solution also provides the 
most consistent overall picture. The spot is similar in 
temperature to the primary {Tspot = 37,560 K), just 
as we would expect if it originates from heating, and 
both the secondary and combined secondary-|-spot tem¬ 
peratures {T 2 +spot = 32, 615 K) agree with the value ob¬ 
tained from the less compromised secondary equivalent 
width ratios. Finally, the combined secondary-|-spot and 
primary temperatures are in better agreement with the 
temperature difference implied by the light curves. This 
solution yields a distance of 1.32 kpc. 

However, if we accept the temperatures estimated 
strictly from equivalent widths, (i.e., Ti = 34,500 K and 
T 2 = 32,100 K), then our best fit is a double contact con¬ 
figuration with a spot inclined ~23° from the northern 
pole, centered ~20° in the direction of rotation from the 
point of contact, and with a radius of ^49°. The average 
secondary-l-spot temperature is ^34,800 K. This does not 
agree with the significant difference in component tem¬ 
peratures suggested by the light curve. Additionally, the 
spot temperature required with this geometry is nearly 
45,000 K, which may be possible with a mass stream 
origin but not heating. The computed distance for this 
scenario is 1.27 kpc. 

The computed distance to all explored configurations 
was between 1.21 kpc and 1.36 kpc, with the closest dis¬ 
tance resulting from a solution with Ti = 34, 000 K, 


T 2 = 28, 596 K, and Tspot = 34,315 K. In all solutions, 
most light-curve-derived orbital parameters had little 
variation (e.g., i ~ 59°, ~ 16.5 R©, and R 2 — 11 R©)- 

The largest source of variation/uncertainty originates 
from the temperatures and resultant bolometric mag¬ 
nitudes, with the second largest source being the high 
uncertainty in interstellar reddening. Because of the de¬ 
generacies present in this system, when computing un¬ 
certainties, we folded in the full variation of all values 
for all explored configurations. We provide these with 
the full list of orbital parameters in Table 2 and adopt a 
distance of 1.32±0.13 kpc to Schulte 3. Additionally, the 
light curve residuals are systematic beyond the stated un¬ 
certainties, indicating additional physical effects that are 
not correctly reproduced in any of the explored models. 

5. DISCUSSION 

The weighted mean distance of the four eclipsing bi¬ 
naries analyzed here is 1.33±0.06 kpc. The unweighted 
mean distance is 1.35 kpc. This value is broadly consis¬ 
tent with the 1.3 - 1.5 kpc distances measured for radio 
masers in pre-main-sequence objects in the surrounding 
molecular clouds (Rygl et al. 2012). It is also in good 
agreement with the Dzib et al. (2012) estimate of 1.3 - 
1.4 kpc toward the massive binary Schulte 5. However, 
our mean distance is somewhat smaller than many mea¬ 
surements in the literature which lie in the range 1.4 - 
1.7 kpc. In particular, Hanson (2003) obtained a spec- 
trophotometric distance of 1.45 kpc toward an ensemble 
of Cyg OB2 stars. Interestingly, Hanson (2003) noted 
that had she adopted the (at that time, newly revised) 
cooler stellar effective temperature scale of Martins et 
al. (2002, 2005) this would place Cyg OB2 at 1.2 kpc 
(D.M.=10.44). She rejected this possibility, concluding 
that such a small distance would reduce the luminosity of 
Cyg OB2 O supergiants that have known mass loss rates 
and create a problem for the theoretical understanding 
of the mechanisms that drive massive star winds. How¬ 
ever, given recent downward revisions of the stellar mass 
loss rates (Puls et al. 2008) it is no longer clear that this 
is a compelling reason for rejecting a smaller distance. If 
we were to adopt the Hanson (2003) distance of 1.20 kpc 
and correct it upward by a factor of 1.08 to account for 
the additional luminosity of secondary stars in close mas¬ 
sive systems as recommended by Kiminki & Kobulnicky 
(2012), the resulting distance is 1.30 kpc, in excellent 
agreement with our current eclipsing binary estimate. 

We note here that Wright et al. (2015) derived 
slightly lower extinctions of Ay = 6 . 98 ^ q ;26 
5.8llg;29 mag for MT91 372 and MT91 696, respectively, 
by using a different set of UV photometry. If we had 
adopted these extinctions, it would have increased the 
resulting distances to these systems to 1.39 and 1.34 kpc, 
respectively. The overall weighted mean distance would 
have risen slightly to 1.36 ±0.10 kpc. This small system¬ 
atic difference is within the uncertainty of our measure¬ 
ment, but it underscores the sensitivity of the distance 
measurement to the adopted correction for interstellar 
extinction. 

The level of agreement between our four distance mea¬ 
surements is remarkable given the stated (perhaps overly 
generous) uncertainties stemming primarily from uncer¬ 
tainties on effective temperatures and interstellar redden- 
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ings. We regard this to be a coincidence given that the 
analysis of each system has entailed independent mea¬ 
surements of these key observable quantities arising from 
independent datasets. The adopted photometry, redden¬ 
ings, and stellar effective temperatures are independently 
justified in each case, and we do not believe that an ar¬ 
tificial level of agreement has been imposed during the 
analysis process. 

Figure 13 displays a three-color representation of the 
Cyg OB2 vicinity, with blue, green, red depicting the 
Spitzer Space Telescope IRAC 4.5 ^m, IRAC 8.0 pm, 
and MIPS 24 pm bands, respectively. Figure 13 shows 
the complex nature of this region, including stellar pho¬ 
tospheres (blue), photo-dissociation regions at the edges 
of molecular clouds (diffuse green tracing broad emission 
features arising from excited polycyclic aromatic hydro¬ 
carbons, and hot dust (red). White circles indicate the 
48 known binary systems in Cyg OB2 (Kobulnicky et al. 
2014), while magenta labels highlight the known eclipsing 
binary systems and their distances. The yellow scale bar 
indicates the linear extent of 10 pc at 1.33 kpc distance. 
Blue labels mark four masers having radio parallax dis¬ 
tances (Rygl et al. 2012) and the massive eclipsing binary 
Schulte 5 (Dzib et al. 2012). 

The maser distance measurements range from 1.3 - 1.5 
kpc, in agreement with most of the eclipsing binaries. 
Given the large angular extent of the Cygnus X region 
(almost a degree!), the entire star-forming complex con¬ 
taining Cyg OB2 could reasonably stretch 100 pc or more 
along the line of sight, thereby encompassing the various 
massive young stellar objects with maser parallaxes. The 
eclipsing binary CPR2002 B17 has a rather uncertain 
distance at 1.5 - 1.8 kpc (Stroud et al. 2010) and could 
be part of an extended star-forming complex somewhat 
further away than Cyg OB2. Schulte 5 is something of 
an enigma. Both Linder et al. (2009) and Ya§arsoy & 
Yakut (2014) analyze the joint light-curve and velocity- 
curve solution to compute distances of 0.93 and 0.97 kpc, 
respectively. While they both present evidence for lumi¬ 
nosity from a third body in this system, the magnitude of 
this contribution is not known. Nevertheless, a third or 
fourth body in the Schulte 5 system (Linder et al. 2009; 
Kennedy et al. 2010) would have to contribute nearly as 
much light as the two eclipsing components in order to 
raise the distance to the mean of other Cyg OB2 eclipsing 
systems. The radio parallax measurement of Schulte 5 
by Dzib et al. (2012) yield 1.3 - 1.4 kpc, in good agree¬ 
ment with the other eclipsing binaries. Given that there 
is no known population of massive stars substantially in 
the foreground to Cyg OB2, we think it unlikely that 
Schulte 5 is a foreground object. Neither is it likely to 
be a runaway in the radial direction; its systemic velocity, 
reported as -55 km s“^ (Ya§arsoy & Yakut 2014), would 
allow it to travel <200 pc in the ^3 - 4 Myr lifetime of 
the most massive stars. A refined understanding of the 
components of Schulte 5 are likely to yield a larger dis¬ 
tance in better agreement with the masers and eclipsing 
systems. 

6 . CONCLUSIONS 


We have presented a photometric and spectroscopic 
analysis of four double-lined eclipsing binaries within the 
Cygnus OB2 Association. The joint analysis provides a 
measure of the temperatures, luminosities, and thereby, 
the distances to these systems. We find distances of 
1.33±0.17, 1.32±0.07, 1.44±0.18, and 1.32±0.13 kpc 
toward MT91 372, MT91 696, CPR2002 A36, and 
Schulte 3, respectively. We adopt a weighted mean of 
1.33±0.06 kpc as the best estimate for the distance to 
the Cygnus OB2 Association. These are the most direct 
distance estimates to this touchstone region of massive 
star formation, and they agree well with radio VLBI par¬ 
allax measurements toward masers in the star-forming 
clouds on the periphery of Cyg OB2. Having a secure 
distance toward the stellar population of Cyg OB2 al¬ 
lows more certainty in the luminosities, mass loss rates, 
and other fundamental stellar quantities measured for 
the massive stars therein. Within several years we an¬ 
ticipate parallaxes from the Gaia space mission provid¬ 
ing even more precise distances toward Cygnus X con¬ 
stituents. Although Gaia’s precision is likely to be lim¬ 
ited in crowded regions near the Galactic Plane, refined 
post-mission analysis of the data may yield a parallax 
distance to a prominent region such as Cygnus OB2. It 
is our hope that this eclipsing binary distance may serve 
as a benchmark for such an attempt. 
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APPENDIX 

TABLE A.l 

Radial Velocities for MT91 372 & MT91 696 


Date (HJD-2,400,000) 

0 

Vrl 

(km s~^) 

Oi 

(km 

-Cl 

s-1) 

Pr2 

(km s“^) 

O 2 -C 2 
(km s~^) 

MT91 372 

54,754.643. 

0.633 

79.8 (14.2) 


-16.5 

-258.2 (14.4) 

-19.0 

54,756.629. 

0.124 

19.4 (17.2) 


3.1 



54,758.616. 

0.416 

-72.0 (21.8) 


5.2 

181.5 (36.2) 

28.5 

56,437.744. 

0.237 

-139.2 (8.9) 


-2.4 

295.8 (11.5) 

5.9 

56,438.831. 

0.725 

96.8 (10.0) 


-26.8 

-278.5 (17.6) 

23.8 

56,439.755. 

0.140 

-120.1 (10.2) 


-17.9 

221.6 (13.8) 

10.0 

56,440.777. 

0.199 

71.1 (9.9) 


-4.5 

-177.4 (12.3) 

14.8 

56,445.838. 

0.871 

77.9 (5.7) 


-5.3 

-208.1 (9.5) 

2.9 

56,446.752. 

0.281 

-150.7 (9.0) 


-12.8 

306.0 (12.4) 

13.8 

56,458.825. 

0.701 

127.0 (12.7) 


6.4 

-269.4 (15.6) 

26.0 

56,459.828. 

0.151 

-122.3 (8.2) 


-14.1 

242.6 (11.6) 

17.5 

56,466.865. 

0.311 

-137.2 (8.1) 


-4.6 

273.7 (12.5) 

-6.1 

MT91 696 

54,699.741. 

0.288 

-234.4 (9.2) 


7.0 

302.9 (14.8) 

-0.4 

54,700.772. 

0.990 

-4.9 (12.3) 


-22.2 

11.5 (15.5) 

30.6 

54,700.902. 

0.078 

-139.4 (22.9) 


-23.0 

150.1 (27.6) 

2.7 

54,701.703. 

0.623 

187.9 (7.3) 


12.4 

-233.9 (13.7) 

-17.7 

54,701.910. 

0.764 

249.1 (9.4) 


-0.4 

-319.9 (13.6) 

-11.3 

54,724.617. 

0.220 

-247.4 (5.1) 


-3.5 

300.0 (8.6) 

-6.4 

54,724.728. 

0.295 

-235.5 (8.8) 


2.7 

296.5 (10.9) 

-2.8 

54,724.824. 

0.360 

-178.1 (19.4) 


12.6 

225.5 (29.0) 

-14.5 

54,725.609. 

0.895 

129.0 (6.1) 


-25.2 

-166.7 (8.3) 

23.0 

54,725.736. 

0.981 

3.5 (8.8) 


-26.9 

17.3 (10.6) 

52.8 

54,725.854. 

0.061 

-114.5 (7.8) 


-21.7 

122.5 (9.8) 

4.4 

54,726.629. 

0.189 

137.4 (6.5) 


4.1 

-181.4 (11.9) 

-17.8 

54,726.743. 

0.667 

218.0 (7.3) 


0.4 

-282.9 (9.8) 

-14.1 

54,726.842. 

0.734 

255.4 (16.8) 


6.3 

-318.7 (21.9) 

-10.6 

54,729.629. 

0.631 

190.0 (6.6) 


6.2 

-257.7 (10.1) 

-31.1 

54,729.724. 

0.696 

242.2 (5.9) 


6.2 

-302.2 (7.8) 

-10.1 

54,729.813. 

0.757 

255.0 (6.9) 


4.8 

-297.6 (8.8) 

11.7 

56,218.613. 

0.111 

-162.4 (8.9) 


-3.0 

200.1 (11.5) 

-0.9 
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Vrl 

Oi - Cl 

Er2 

02 -C2 

Date (HJD-2,400,000) 

0 

(km s“^) 

(km s“^) 

(km s“^) 

(km s“^) 

56,222.568. 

0.803 

231.6 (5.9) 

-5.1 

-280.3 (8.2) 

12.2 

56,218.632. 

0.124 

-167.3 (10.9) 

7.2 

242.6 (13.7) 

22.8 

56,222.788. 

0.953 

62.6 (7.9) 

-11.1 

-85.1 (11.8) 

4.3 

CPR2002 A36 

54,403.610. 

0.162 

-153.8 (7.1) 

5.3 

177.8 (4.2) 

1.6 

54,403.702. 

0.181 

-167.8 (6.7) 

-2.3 

193.8 (3.9) 

2.8 

54,405.705. 

0.610 

9.9 (15.4) 

-55.3 

-208.5 (9.1) 

-20.3 

54,406.632. 

0.808 

110.2 (4.8) 

-0.1 

-257.5 (2.9) 

-9.0 

54,406.757. 

0.835 

104.7 (6.2) 

2.2 

-243.7 (3.7) 

-10.7 

54,408.664. 

0.243 

-185.7 (11.0) 

-12.0 

234.9 (6.5) 

20.7 

54,409.730. 

0.471 

-39.9 (48.9) 

16.0 

-14.5 (28.8) 

-50.6 

54,410.667. 

0.671 

90.8 (7.1) 

-7.9 

-245.0 (4.2) 

-7.8 

54,641.827. 

0.119 

-100.3 (10.2) 

38.4 

121.2 (6.0) 

-11.2 

54,642.754. 

0.317 

-176.2 (7.1) 

-15.7 

211.5 (4.2) 

16.8 

54,644.770. 

0.748 

116.9 (9.1) 

1.0 

-252.1 (5.4) 

10.3 

54,645.834. 

0.976 



-25.6 (6.0) 

47.4 

54,647.767. 

0.389 

-148.2 (11.5) 

-25.1 

115.4 (6.8) 

-19.8 

54,669.843. 

0.112 

-103.4 (8.7) 

31.7 

121.5 (5.1) 

-2.0 

54,672.777. 

0.739 

109.5 (7.2) 

-5.6 

-276.0 (4.2) 

-14.3 

54,673.875. 

0.974 


-62.3 (27.4) 

13.3 

54,696.684. 

0.853 

121.3 (6.8) 

26.4 

-229.1 (4.0) 

-9.7 

54,696.964. 

0.913 

101.5 (8.0) 

37.1 

-158.7 (4.7) 

-0.6 

54,697.659. 

0.062 

-68.8 (19.5) 

37.4 

14.0 (11.5) 

-41.6 

54,698.753. 

0.296 

-196.2 (11.4) 

-29.5 

206.3 (6.7) 

1.0 

54,700.871. 

0.749 

103.6 (9.0) 

-12.4 

-259.7 (5.3) 

2.7 

54,701.895. 

0.968 


-78.4 (29.5) 

6.3 

Schulte 3 

53,989.652. 

0.643 

42.1 (5.7) 

15.6 

-224.4 (3.4) 

0.4 

53,989.774. 

0.668 

37.4 (4.1) 

3.4 

-238.6 (3.4) 

5.8 

53,990.848. 

0.895 

16.3 (7.2) 

-0.8 

-191.9 (4.0) 

-3.5 

54,286.904. 

0.276 

-145.4 (8.3) 

-7.2 

201.2 (8.4) 

28.0 

54,341.802. 

0.843 

117.2 (12.3)“^ 

82.6 

-217.4 (5.4) 

19.1 

54,342.724. 

0.038 

5.0 (104.8E 

82.0 

13.1 (46.0) 

8.1 

54,343.830. 

0.271 

-134.0 (6.0E 

4.7 

249.3 (2.6)“" 

75.2 

54,344.761. 

0.467 

-21.1 (36.4E 

46.2 

11.0 (16.0) 

1.0 

54,345.781. 

0.682 

90.1 (9.8E 

53.0 

-251.2 (4.3) 

1.5 

54,346.787. 

0.894 

54.5 (25.4)“* 

37.0 

-170.1 (11.1) 

19.3 

54,347.755. 

0.098 

-19.0 (9.7E 

89.8 

101.4 (4.3) 

19.3 

54,348.840. 

0.326 

-116.8 (5.7E 

13.1 

220.4 (2.5)“" 

66.8 

54,628.748. 

0.305 

-128.3 (4.7) 

6.1 

172.4 (5.8) 

8.0 

54,628.880. 

0.333 

-134.2 (4.0) 

-5.8 

162.7 (3.3) 

13.1 

54,629.866. 

0.141 

-49.7 (3.3) 

-26.8 

-126.7 (3.2) 

-16.7 

54,630.739. 

0.725 

32.5 (4.4) 

-11.1 

-260.6 (2.6) 

8.7 

54,630.852. 

0.749 

40.2 (4.1) 

-4.8 

-258.6 (2.8) 

13.5 

54,632.745. 

0.147 

-122.1 (4.5) 

4.7 

155.4 (4.5) 

24.2 

54,633.741. 

0.357 

-140.0 (5.6) 

-18.9 

128.0 (5.9) 

-3.7 

54,633.859. 

0.382 

-120.0 (3.0) 

-8.3 

82.9 (3.6) 

-26.7 

54,642.717. 

0.249 

-128.9 (6.4)“" 

10.9 

258.2 (2.8)“" 

82.9 

54,642.878. 

0.283 

-143.6 (3.8)“" 

-6.1 

252.9 (1.7)“" 

81.0 

54,644.684. 

0.663 

91.2 (6.4)“" 

58.6 

-254.0 (2.8) 

-13.3 

54,645.734. 

0.884 

101.5 (13.8)“" 

80.3 

-188.2 (6.1) 

11.6 

54,646.707. 

0.089 

-39.9 (8.9)“" 

65.3 

100.8 (3.9) 

28.7 

54,647.688. 

0.296 

-139.1 (8.0)“" 

-3.2 

229.0 (3.5)“" 

61.2 

54,648.902. 

0.152 

-13.9 (12.1)“" 

2.6 

-115.8 (5.3) 

10.0 

54,672.692. 

0.165 

0.3 (5.9)“" 

9.6 

-142.7 (2.6) 

-0.1 

56,218.402. 

0.257 

-124.8 (8.0)“" 

14.6 

223.8 (3.5)“" 

47.4 

56,219.341. 

0.455 

-56.1 (51.0)“" 

17.5 

-22.7 (22.4) 

-48.4 

56,222.349. 

0.089 

-81.9 (7.9)“" 

24.0 

172.7 (3.5)“" 

98.1 


Observation not used in solution fit. 
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TABLE A.2 

V-BAND PHOTOMETRY OF MT91 372 


Kiminki et al. 


TABLE A.4 

V-BAND PHOTOMETRY OF CPR2002 A36 


HJD-2,400,000 

V 

(days) 

(mag) (mag) 


56,458.819. 

. 14.983 

0.011 

56,458.821. 

. 14.991 

0.012 

56,458.822. 

. 14.936 

0.011 

56^458.825. 

. 14.955 

0.011 

56^458.827. 

. 14.962 

0.011 

56,458.830. 

. 14.977 

0.011 

56^458.832. 

. 14.993 

0.011 

56^458.837. 

. 14.986 

0.011 

56,458.839. 

. 14.994 

0.011 

56(458.841. 

. 14.964 

0.011 


HJD-2,400,000 

V 

(days) 

(mag) (mag) 


56,271.541. 

. 11.421 

0.002 

56,271.542. 

. 11.432 

0.002 

56,271.543. 

. 11.428 

0.002 

56,271.648. 

. 11.447 

0.002 

56,271.649. 

. 11.434 

0.002 

56,271.649. 

. 11.429 

0.002 

56,271.650. 

. 11.435 

0.002 

56,271.650. 

. 11.428 

0.002 

56,271.651. 

. 11.430 

0.002 

56,271.651. 

. 11.439 

0.002 


Note. — Table A.2 is published in its entirety in the electronic 
edition of the Astrophysical Journal. A portion is shown here for 
guidance regarding its form and content. 

^ 1 (7 error includes formal photon statistics uncertainties only; addi¬ 
tional uncertainties at the level of 0.01 mag are present as a result of 
uncorrected atmospheric transparency variations. 


TABLE A.3 

V-BAND PHOTOMETRY OF MT91 696 


HJD-2,400,000 

(days) 

V 

(mag) 

(T® 

(mag) 

56,219.567. 

... 12.066 

0.002 

56,219.569. 

... 12.068 

0.002 

56,219.570. 

... 12.064 

0.002 

56,219.572. 

... 12.067 

0.002 

56,219.573. 

... 12.067 

0.002 

56,219.589. 

... 12.065 

0.002 

56,219.590. 

... 12.069 

0.002 

56,219.592. 

... 12.063 

0.002 

56,219.593. 

... 12.067 

0.002 

56,219.595. 

... 12.069 

0.002 


Note. — MT91 696 photometry includes the contribution from the 
blended 0'/02 tertiary MT91 696c discussed in the text. Table A.3 is 
published in its entirety in the electronic edition of the Astrophysical 
Journal. A portion is shown here for guidance regarding its form and 
content. 

^ 1 (7 error includes formal photon statistics uncertainties only; addi¬ 
tional uncertainties at the level of 0.01 mag are present as a result of 
uncorrected atmospheric transparency variations. 


Note. — Table A.4 is published in its entirety in the electronic 
edition of the Astrophysical Journal. A portion is shown here for 
guidance regarding its form and content. 

^ 1 <7 error includes formal photon statistics uncertainties only; addi¬ 
tional uncertainties at the level of 0.01 mag are present as a result of 
uncorrected atmospheric transparency variations. 

TABLE A.5 

V-BAND PHOTOMETRY OF SCHULTE 3 


HJD-2,400,000 

(days) 

V 

(mag) 

(mag) 

56,218.586. 

... 10.194 

0.001 

56,218.587. 

... 10.182 

0.001 

56,218.587. 

... 10.169 

0.001 

56,218.588. 

... 10.168 

0.001 

56,218.588. 

... 10.171 

0.001 

56,218.589. 

... 10.175 

0.001 

56,218.589. 

... 10.181 

0.001 

56,218.590. 

... 10.156 

0.001 

56,218.590. 

... 10.174 

0.001 

56,218.591. 

... 10.149 

0.001 


Note. — Schulte 3 photometry includes the contribution from the 
blended tertiary Schulte 3C discussed in the text. Table A.5 is pub¬ 
lished in its entirety in the electronic edition of the Astrophysical Jour¬ 
nal. A portion is shown here for guidance regarding its form and con¬ 
tent. 

^ 1 (7 error includes formal photon statistics uncertainties only; addi¬ 
tional uncertainties at the level of 0.01 mag are present as a result of 
uncorrected atmospheric transparency variations. 




























































